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1 Introduction 

Identical particles with statistics continuously interpolating between Bose-Einstein and 
Fermi-Dirac statistics can exist in one and two dimensions |l|]. In one dimension, these 
statistics are described by the Calogero model , in the sense that if the particles are not 
classically free, their asymptotic properties are however not affected by the interaction, 
up to a permutation. In two dimensions, they are described by the anyon model: in 
the singular gauge, free anyonic eigenstates pick up a phase vra when any pair of anyons 
are exchanged [Q, where a is the statistical parameter. The spectrum of the Calogero 
model in a harmonic well is quite similar to those of the anyon model in a harmonic well 
projected on the lowest Landau level (LLL) of an external magnetic field [^. 

It is now widely accepted that particles with intermediate (fractional) statistics should 
play a role in the fractional quantum Hall effect 0, a strong magnetic field, low temper- 
ature effect observed in certain bidimensional conductors. When an integer fraction 1/m 
of the LLL is filled by the electrons, the Hall conductivity exhibits a plateau at a value 
Xjm in unit of hc/e^, and the longitudinal resistivity vanishes. The trial eigenstates pro- 
posed by Laughlin for the groundstate and the excited states are sustained by numerical 
few-body computations @]. On the Hall plateaux, the Coulomb repulsion lifts the Lan- 
dau degeneracy and leads to a nondegenerate groundstate with a gap, which explains 
the absence of dissipation in the longitudinal transport. The excitations are described 
as quasiparticles or quasiholes localised on some defects of the sample. A Berry phase 
calculation shows that they obey intermediate statistics, i.e. they are anyons. More re- 
cently, these excitations have been shown to obey a generalisation of the Pauli principle 
first proposed by Haldane ||^, ^J. The general relation between the anyon model and this 
generalised exclusion principle remains to be clarified. It has however been shown that it 
is obeyed by anyons in a strong magnetic field at low temperature [^, |TU . 



On the one hand, the statistical mechanics of an anyon gas in a strong magnetic field at 
low temperature has been analyzed |]1T| in a situation where the singular flux tubes carried 



by the anyons are anti parallel to the magnetic field. In this particular screening regime. 



the Hilbert space is restricted to the Landau groundstate, with effectively one quantum 
number per particle, its orbital angular momentum. There exists a critical value of the 
filling factor, where the screening is complete, and for which the system is incompressible 
(nondegenerate with a gap), leaving a Bose condensate. 

On the other hand, the statistical mechanics of the Calogero model has been recently 



analyzed by various ways |jT2[, [T^. 



The present letter presents a thorough discussion of the anyon model in a strong 
magnetic field, and displays its intimate relation with the other concepts of intermediate 
statistics. 

In section 2, the cluster coefficients for the Calogero and anyon models will be com- 
puted in a harmonic well, and will be shown to have the same form. Here, the harmonic 
well has to be understood as a long distance regulator, which, when it is properly taken 
to vanish, yields the correct thermodynamic limit. Indeed, in sections 3 and 4, the ther- 
modynamic limit of the anyon and Calogero models will be considered. The general 
thermodynamical prescription, which is depending on the space dimension d, is that, at 
order n in the power series expansion of the thermodynamical potential in the fugacity 
2, one replace l/^njS'^uj'^Y^'^ — > V/X^ {V is the volume of the system and A = J2Txl3/m 
is the thermal wavelength) |T^. An appendix follows where the thermodynamic limit of 
the multi- species statistical problem is considered. 

2 Statistical mechanics of the anyon and Calogero 
models in a harmonic well 

In the sequel, one will consider Bose-based statistics quantum mechanics : eigenstates 
will be symmetric and the statistical parameter a by convention such that a = for 
Bose statistics, and a = —1 for Fermi statistics. In the anyon model, a measures the 
statistical flux carried by the anyons in unit of the flux quantum. The system is periodic 
in a with period 2. In the Calogero model, a is related to the coupling constant g of the 



one dimensional Calogero interaction g/xfj hy g = a{a + 1). 

In a harmonic well tu, the spectrum of the anyon model in the LLL of an external 
magnetic field and the spectrum of the Calogero model are similar. As for ideal bosons, 
an A^-body eigenstate happens to be entirely characterized by the number n^ of 1-body 
eigenstates with a given orbital quantum number £ = 0, 1, ...oo and energy e^ = eo + ^ro, 
with the constraint ^i rig = N. The A^-body spectrum is nothing else but the sum of the 
1-body spectrum X^^n^e^ shifted by the 2-body interaction energy —^N{N — l)rz!a. The 
anyon spectrum in the LLL corresponds to tu = ujt—ujc (i^t = J^c + ^"^^ ^c = \^B\/{2m)) 
and eo = LJt, whereas the Calogero spectrum corresponds to cu = cu and eo = juj. 

Thus, the net effect of a intermediate statistics is simply the shift of the A^-body Bose 
spectrum by —^N{N — l)zua, which in turn affects the Boltzmann weight of the A^-body 
bosonic partition function 

Zn = e5/3^(^-i)-" TT -^ (1) 

The cluster expansion of the thermodynamical potential fi = — IuX^at Zj^jz^ = — X^n ^n-2" 
in power series of the fugacity z = exp /?/i, where /i is the chemical potential, yields, in 
the limit where the harmonic attraction becomes small, i.e. vu ^ 0, ITTI, |T2ll 



^^^^^^^n^ (2) 

Clearly, the cases a = and a = —1 coincide with the bosonic and fermionic cluster 
expansions. Note also that the polynomial expression (j^) appears in some recent numerical 
estimation of the A^-cycle brownian closed paths contribution to the A^-anyon partition 
function |]T3 . 



Let us consider the probability V{nf) to have a given occupation number ni . For 
bosons in a harmonic well, it is well known that Vh{ne) = (1 — ze~^^''){ze~^'^'')"''^. For 
particles with a statistics in a harmonic well, it should be defined as 

PK)=e^tr„,^e-^(^--'^^) (3) 

where the trace is as usual understood as summation over the n/s but with the constraint 
that one of them, rig^, is fixed. 



The trace can be computed from its bosonic end value e~^^Vh{ne) by the shift z'^ —>■ 
^^f3N{N-i)u^a^N ^^ ^^^q power series expansion in z. In terms of the intermediate statistics 
grand partition function e~^^^\ one deduces 

In the limit xu ^ 0, one approximates 

oo 

- lim (fi(ze"^^^") - n{z)) = ^ n^aPwnbnZ'' = n^aXny (5) 

n=0 

where y is solution of 

y - ze-^'^'y^^'' = 1 (6) 

with y — > 1 when z -^ []TB[. One directly finds 



1 



The occupation numbers are essentially statistical independent in the approximation used 
in deriving (7). Indeed, the joint probability V{ni,nm) ~ V{ni)V{nm) factorizes, which 
implies subleading correlations {neUm) — {ni) (nm) compared with {ne){nm)- Alternatively, 
correlations become important when a large number of occupation numbers is considered. 

Note that one could as well have dealt with Fermi-based statistics quantum mechanics, 
meaning antisymmetric eigenstates and fermionic occupation numbers ne either or 1. 
One then would have found V{ne) = (y-l)"^e-"^'^('*"'"V(l + {y - l)e-^('*~'o)) and {ne) = 
V{1). 

One is now in position to consider the thermodynamic limit cj — > in (1,2). The 
prescription has already been given in the introduction. It is simply dictated by the 
natural requirement that when u ^ 0, the thermodynamical potential of a system of 
(i- dimensional harmonic oscillators coincide with those of a system of d-dimensional par- 
ticles in a box of infinite volume. A qualitative justification is that in a neighbourhood of 
any given point r, one can always approximate X^i \rnuj'^r'f by ymtu^r^. It follows that for 



a system of particles confined in an infinitesimal volume df'r around r, the local thermo- 
dynamical potential in the presence of a small harmonic attraction can be approximated 
by simply replacing z —>■ ze'^^"^^ "^ and V —>■ (Pr in the infinite box thermodynamical 
potential. Since it is additive in the limit tu — * 0, for the entire system one has 

— Tb(z e-|/3--'-'r = - V - b z^ (9) 

where d is the space dimension of the system. Thus the prescription l/(n/5^ci;^)'^/^ — > V/X^ 



at each order 2;" [|1J], that will be used in the following sections. 



3 Statistical Mechanics of the Anyon gas in a strong 
i^-field : the thermodynamic hmit 

The model : The A^-anyon Hamiltonian {h = 1 and c = 1) in an external magnetic 

field B is 

^ ^1/19 ^kxn, B-> \^ 

HN = 2^7^[TT7^~a}_^ ^-e—kxn) 10 

fr[2m\idri ^. rfj 2 / 

where k is the unit vector perpendicular to the plane, and fij = fi—fj. The Hamiltonian 

being invariant under {xi,yi,a,e) — ^ (xj, — j/j, —a, — e), where e = eB/\eB\, the spectrum 

and thus the partition function are invariant under (a, e) -^ (— «, — e). They depend only 

on |a|,ea and uJc- By convention, and without loss of generality, one chooses e = +1 

(in the opposite case one would simply change a —>■ —a). The shift a —>■ a + 2 in H]\f 

amounts to the regular gauge transformation ip -^ exp{—2iJ2i<j argrjj)?/', which does not 

affect the symmetry of the eigenstates, implying that the spectrum is periodic in a with 

period 2. 

One finds an infinitely degenerate groundstate with energy Nuc and eigenstate basis 

(in the interval a G [—1,0]) |]lj 



^ = Er-^'"^E4' eM-lniuj,J2z,z,), e, > o (ii) 



i<j i 



where 5 is a symmetrisation operator. If one leaves aside the anyonic prefactor Y[i<j '^ij°'^ 
the iV-anyon groundstate is a symmetrised product of 1-body Landau groundstates of 



energy ujc and orbital angular momentum £j. One can argue that, starting from the Bose 
case a = modulo 2, the gap above the groundstate remains of order the cyclotron gap 2uJc 
when a G [—1, 0] modulo 2, but is no more under control in the interval a G [0, 1] modulo 
2. In the limit a — > 0"^ modulo 2, unknown non-linear states (only a small subset of excited 
states are known) have to merge in the groundstate basis to get a complete Landau basis. 
This result is sustained by various numerical and semi-classical considerations. It could 
even be true, as a numerical computation of the 4-anyon spectrum in a magnetic field 
seems to indicate, that the gap above the groundstate is exactely equal to the cyclotron 



gap for a G [—1,0] modulo 2 [|18|. Thus when a G [—1,0] interpolates between Bose 
and Fermi statistics with singular flux tubes anti-parallel to the external 5-field, the 
thermal probability e"^^'^'' to have an excited state is negligible when the thermal energy 
kT = 1/(3 is much smaller than the cyclotron gap. In what follows, one focuses on the 
thermodynamics of an anyon gaz precisely in these conditions where it is licit to consider 
anyons projected in the groundstate of the external magnetic field. 

Statistical mechanics in the thermodynamic limit : In order to define a proper 
thermodynamic limit, i.e. to give a non ambiguous meaning to the infinite degeneracy of 
the spectrum and to the factorized volume in the thermodynamical potential, the system 
should be regularized at long distance. With a harmonic regulator, the eigenstates are still 
given by (0), but with c^c — * ^t- The harmonic confinement partially lifts the degeneracy, 
and, as already emphasized, the intermediate a statistics shifts the A^-body Bose spectrum 
J2i Q, by -|A^(iV - l){uJt - ujc)a, where the 1-body energy is e^^ = ujt + iiiujt - cUc). In 
2 dimensions, the thermodynamic limit uj -^ in the cluster expension is l/{n(3^uj^) — > 
V/X^. One precisely finds a thermodynamical potential [11 

n = -pi^Vlny (12) 

where pl = 2/5ciJc/A^ is the Landau degeneracy per unit volume, and y is the implicit 
solution of equation (|^) with eo = uJc- One can check that the ideal bosonic and fermionic 
thermodynamical potentials are readily recovered when a = and a = — 1, since y = 



1/{1 — ze ^"^t:), and respectively y = 1 + ze ^^''. The ratio of the mean anyon number ^ 
to the Landau degeneracy is given by 

A^L ^^ -a l + av 

u is monotically increasing with z, as required. Its imphcit definition directly follows 

(1 + i^ + az/)i+"(l + az/)-° ^ ' 

Standard thermodynamical functions as the pressure P = —Q/{Vl3), the magne- 
tization Ai = —dQ/dB, the internal energy U = dQ/df3 + fi{N) and the entropy 
TS = U — fJ,{N) — Q/P follow from (|12j) and (|T3|) : the equation of state reads 



and the virial coefficients are 



P/3 = pLln(l + -^^) (15) 

1 + auy 



^ ^"-^^-{(1 + «)"-«"} (16) 



PL n 

In the limit where the Boltzman weight exp{—2j3ujc) can be neglected, both the exact 
second virial coefficient and the pressure at the first perturbative order in a |TP 



are 



correctly reproduced. We will comment later on the divergence of the pressure at the 
critical value u^ = — 1/a. 

The magnetization per unit volume is 

M = -PBP + 2^ In (l + -^) (17) 

where fi-Q = \e\/2in is the Bohr magneton. Except near the singularity v = —1/a, the 
ratio of the logarithmic correction to the De Haas- Van Alphen magnetisation —^bP is of 
order (/JtUc)"^, and thus negligible. 

The internal energy U = {N)uJc and the entropy S reads 

TS = ^piy In ^ "^^ [,^ (18) 



^It can be verified a posteriori that the fluctuations {{N'^) — (N)'^) / {N^ vanish in the thermodynamic 
limit. 
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The expression pq ) for the entropy has a natural interpretation. If one assumes that 
the small 1-body angular momenta £« are energetically preferred and that they are the 
only one to survive in the thermodynamic limit, one infers that the number G of quantum 
states available per particle should be defined by < £j < G. Since the bosonic A^-anyon 
groundstate is the symmetrised product of 1-body eigenstates of momentum £i (|ll]), its 



degeneracy is Cq^^^j^. Accordingly, since the entropy (0) is the logarithmic measure 



of the groundstate degeneracy, one has S = klnCir , q, wAr\, from which one deduces 



{G)c^N^ + a{N-l). 

G also appears in the mean occupation number (rii), simply because it is independent 
on i in the thermodynamic limit. One indeed verifies that (A^) = J2e{''T-e) = {G){n() 
implies (G) = A^l(1 + av). 

Note also that if one assumes that the A^-th anyon experiences an uniform effective 
magnetic field Sefr = B + acpo^N — 1)/V^, i.e. the external magnetic field screened by the 



flux tubes carried by the other A^ — 1 anyons |2y], one simply recovers A^£ — ^ B, 



^cSl<Po 



Ni^ + a{N — 1). In this mean field point of view, the effective magnetic field -Befr vanishes 
at the critical filling. 

The mean field result given above can be strengthened by computing the degeneracy 
per particle directly in a box. To achieve this goal, one determines which eigenstates have 
their maxima contained inside a cicular box of radius R. Let us consider the eigenstates 
(|Tl|) without bothering about symmetrisation to make the argument simpler. At the 
maximum of an eigenstate one has 

^^l^r=f-«E& + 4---crl)|V^r = (19) 

The density of particles is defined as p{r) = {ip\ J^j^i^i^— rj)\tp) / {ip\tp) , normalized to 
/ (Prp{r) = N — 1. It depends only on r because {ipl"^ is rotationally invariant. In a mean 
field approach, one replaces the summation in (p!9[) over discrete indices j by a continuous 
integral over the density 

-a f d^r'p{r ,!'"!,,, + ^i^ ' ^'^-^^ = (20) 

J [Ti — r 'Y rf 



Using the identity di. J d r'p{r '){fi — f ')/(fi — r ') = 27rp(rj), one finds 

ii = mujcrf + a 2nrdrp{r) (21) 

Jo 

One deduces that ii increases with r^, implying a maximum angular momentum corre- 
sponding to the size of the box R. Thus, the maximum of the eigenstate will be inside 
the box if and only if 

rR 

0<ii< mu^R^ + a / 2'Krdrp{r) = A/'l + a{N - 1) (22) 

Jo 

In conclusion, one has 

G = iVL + a{N - 1) (23) 

which correctly reproduces the usual Landau degeneracy when a = (and also when A^ = 
1). The result ( p3[ ) is convincing because it implies that the Bose-based A^-anyon ground- 
state (|TT|) is, in the singular gauge, the Landau groundstate of a fermionic system when 
a = -1, as it should. Indeed, the basis (|TI|) is then {ni<j % S Hi z-'e'^"^" ^» ^*^% < ii < 
Ni^-{N-1)}, whereas the usual Fermi basis is {^ni-zl"'e"5'"'^-E.^«2», < m^ < A^l} {A 
the antisymmetrisation operator) . Both basis are equivalent because any totally antisym- 
metrised polynomial in Zi of degree di can be written as the product of the Jastrow factor 
Y[i<j Zij and of a totally symmetrised polynomial of degree rfj — (A^ — 1), and reciprocally. 

Discussion : (^) coincides with Haldane's statistics definition AG = aAA^ 0, which 
therefore appears as the analytical continuation in a of the anyon statistics in a strong 
magnetic field. Actually, the entropy as well as other thermodynamical quantities can be 
recovered from (19) (see for example |llO|). 

When a = 0, any value of v is allowed, due to Bose condensation. In the case of Fermi 
statistics a = — 1, Pauli exclusion implies that the LLL is completly filled when u = 1. 
At a particular negative a, anyons obey a generalised exclusion principle. Indeed, the 
groundstate exists as long as {G) > 0, i.e. u < —1/a. The critical filling Ucr = — 1/a 
describes a non degenerate groundstate with all the £j's null, i.e. a minimum angular 
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momentum L = — |(A^)cr((^)cr — !)«• In the singular jauge, 

N 

^' = n^.7<^-p(-^E^^^"^) (24) 

i<j i 

when a = —1, one recovers a Vandermonde determinant buih from 1-body Landau eigen- 
states. Incidentally, in the Haldane statistics point of view, when a = —m, the non degen- 
erate groundstate coincides with the Laughlin eigenstate at the critical filling z/cr = l/m. 

Since transitions to excited levels are by construction forbidden, the pressure diverges 
and the entropy vanishes when the LLL is fully occupied (such that any additional particle 
is excluded). The system is in its non degenerate groundstate as one can also see by 
inspecting the large z behavior of the thermodynamical potential f2 — > — e~^^^'*'=' i^c^{N)cT_ 
In this situation the gas is incompressible : the isothermal compressibility coefficient 
Xt = —y (f^j vanishes at the critical filhng. 

On the contrary, when one analytically continuates the thermodynamical quantities 
computed above through positive a, any value of the filling factor is allowed since G 
increases with the filling factor. The pressure becomes constant while the entropy diverges 
when u -^ oo. 

4 Statistical Mechanics of the Calogero model in the 
thermodynamic limit 

As stated in the introduction, the Calogero model can be viewed as an intermediate 
statistics model for 1-dimensional particles. Contrary to the anyon model, there is no 
periodicity in the statistical parameter a. One also assumes a < 1/2 in order to have 
square-integrable eigenstates. The Calogero Hamiltonian is 

A long distance harmonic regularor is again introduced to discretize the spectrum. One 
redefines ip = Y[i<j{xi — Xj)~^^^ where the Hilbert space of eigenstates ip can always be 
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choosen to be completely symmetric. In this way, a new Hamiltonian is obtained [p^ 



H = Y.is4a^ + \)uj - \n{N - l)au (26) 

where the creation and annihilation operators af = -j={\/rriuJXi — Di/^/muJ) and Sj = 
4=(y^mcJxj + Dij ^rauj\ with Di = di — aY^j^ii.^ — Pij)/{xi — Xj) and Pij the exchange 
operator, satisfy the commutation relations [aj,a^] = %(! — aX^^Li -Pifc) + cuPij. The 
groundstate ipo is solution of aiipo = for all i, and the excited states are obtained as 
4'{ei} = 'SYlittj' ^ipo with energy J2ii^i + ^)^ ^ ^N{N — l)auj. Thus, as expected, a 
statistics shifts the Bose spectrum by a constant. 

One would now like to compute the thermodynamical potential in the thermodynamic 
limit, following a procedure analogous to the one used in section 3 for the anyon model. 
However, in one dimension, the thermodynamic limit u ^ has to be understood for 
a given coefficient cluster coefficient bn as l/{y/nPuj) -^ V/\. Using the identity jy^ = 

2 

2^ I^oo ^P e~"'^27^ , one can easily resum the cluster expansion in the thermodynamic limitQ 
to get 

n = -PVP = f dplny (27) 

Here, y is the solution of 

y - ze-^pV^rnyi+a ^ ^ ^28) 

where p^ /2m appears quite naturally as a continuous one dimensional energy spectrum 
in the cu — > limit. The thermodynamical potential for bosons (fermions) is correctly 

2 

reproduced when a = since y = 1/(1 — ze~^2^) (respectively a = —1 since y = 

2 

1 + ze~^2^). One deduces the density 

1 /■ , 1 
p = — dp —. (29) 



27r 



y-i 



a 



and the internal energy 



V r P^ 1 

U = — dp-^ . (30) 

27r y ^ 2m 1 ^ ' 
--a 

y-1 



*ri has also been derived via the Bethe ansatz directly in the thermodynamic limit [|2| p^ . 
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One can also compute the number g{p)dp of quantum states at momentum p available 
per particle and per unit volume. Since the local thermodynamical potential at momentum 
p in (p7D coincides with the anyonic thermodynamical potential ([12D where Uc is replaced 
by p'^/2'm, this suggests that the occupation number Up can be deduced from the anyonic 
one, given in (Q) with u = 0. One infers 

np = y-l (31) 

where y is solution of (pSl). However, ^ J dp Up does not define the total number of 
particles because of possible corrections due to particle correlations. In fact, the informa- 
tion about these corrections is precisely described by g{p). If one reminds that the local 
density of particule at momentum p is nothing but Up g{p)dp, one gets from 



«<P' = £ l-.fa-l) P2) 



which coincides with the "hole density" of the Bethe ansatz \T3]. As expected, the entropy 
can be rewritten as 

S = kl- 1 dp lnC;i^^^^,^, (33) 

Note that the particular thermodynamical limit algorithm in one dimension has played 
a crucial role for the derivation of a thermodynamical potential, which coincides with the 
one obtained directly in the continuum via the Bethe ansatz. Any other prescription to 
define the thermodynamical limit would lead to different results pz 



5 Conclusion and Perspectives 

In conclusion, the Calogero model and the anyon model in the LLL of an external mag- 
netic field have been shown to be intimately connected, not only at the level of their 
spectrum, but also of their thermodynamical properties. By anyon model, one means 
here singular flux tubes anti-parallel to the i?-field, i.e. a G [—1,0]. Haldane's statis- 
tics can be understood in this context, simply as the analytical continuation in a of the 
statistical mechanics of the anyon model in the strong 5-field. 
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It would be certainly interesting to have more precise information on the thermody- 
namics of the anyon model when a G [0, 1]. To do so, one has necesseraly to take into 
account, in one way or an other, the excited states which join the groundstate when a = 0. 
More interesting and richer physics, than the strictly 1-dimensional one, is to be expected 
in this regime. 

Appendix : Many anyons species in the thermodynamic limit 

One can consider different species of identical particles to define a generalized anyon 



model PI 



^ = E ^ ( P^ - E «^.^^ - \^^^^ X ^ (34) 

When all the a^- = a, one recovers the anyon model. One again assumes Cj-B > 0. The 
shift aij — > ajj + l amounts to the regular gauge transformation il) -^ exp(— i J2i<j argry)?/) 
which does not affect the symmetry of the eigenstates. It follows that the spectrum is 
now periodic for each aij with period 1. Using the same lines of reasonning as above, one 
finds the groundstate for a^j g] — 1, 0] {ud = eiBi/2mi) 

^ = n ^T' n (^f exp{-l-muj,iZiZ^)) i, > (35) 

Contrary to anyons, eigenstates are not totally symmetric anymore, implying that some 
excited eigenstates have to merge in the groundstate when a^ — ^ —1. A clear information 
on the gap above the groundstate is missing. However, one can still project the system 
on the LLL Hilbert space. Let us denote the number of particle in the specy a by Aq. 
If one generalizes the mean field computation of the degeneracy per particle given above, 
one recovers the mutual statistics definition of Haldane 

G'a = iVL, + ^a,,(A,-,5,,) (36) 

h 

with Alq = ea-BaV^/27r, where Ga satisfies AG^ = Yl,h'^ah^^h- The total degeneracy is 
thus riaC'c^^.i+jv^ and the thermodynamical potential reads [|I^ 

^=-PVI5 = -Y.N^My^ (37) 
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where Va is solution of i/„ - ze"^'^'^" Tlfo 1/f '^''""'"'^'^'' = 1 with Va ^ I when z ^ 0. The 
filhng factor //„ = {Na)/Ni^a of the specy a is determined by 

ya = l + — ^ (38) 

1 + 2^ aabi^b 



One finds for the magnetization and the entropy 



/"BaPa + 2— r- In 1 + 



A^ V l+Ebttafet'^ 



;i + z/„ + 5]a„5Z/5)i+'^"+S^"-'^'' 



(39) 



S = y kNi^a In ^- -^^ (40) 
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